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Abstract. High resolution, high S/N spectra have been 
obtained for a sample of 90 F and G main-sequence disk 
stars covering the metallicity range —1.0 < [Fe/H] < +0.1, 
and have been analysed in a parallel way to the work of 



Edvardsson et al. (1993a) in order to re-inspect their re- 
sults and to reveal new information on the chemical evo- 
lution of the Galactic disk. 

Compared to Edvardsson et al. the present study in- 
cludes several improvements. Effective temperatures are 
based on the Alonso et al. (1996) calibration of color in- 
dices by the infrared flux method and surface gravities 
are calculated from Hipparcos parallaxes, which also allow 
more accurate ages to be calculated from a comparison of 
My and T e g with isochrones. In addition, more reliable 
kinematical parameters are derived from Hipparcos dis- 
tances and proper motions in combination with accurate 
radial velocities. Finally, a larger spectral coverage, 5600 - 
8800A, makes it possible to improve the abundance accu- 
racy by studying more lines and to discuss several elements 
not included in the work of Edvardsson et al. 

The present paper provides the data and discusses 
some general results of the abundance survey. A group 
of stars in the metallicity range of —1.0 < [Fc/H] < —0.6 
having a small mean Galactocentric distance in the stel- 
lar orbits, R m < 7 kpc, are shown to be older than the 
other disk stars and probably belong to the thick disk. 
Excluding these stars, a slight decreasing trend of [Fe/H] 
with increasing R m and age is found, but a large scatter 
in [Fe/H] (up to 0.5 dex) is present at a given age and R m . 
Abundance ratios with respect to Fe show, on the other 
hand, no significant scatter at a given [Fe/H]. The derived 
trends of O, Mg, Si, Ca, Ti, Ni and Ba as a function of 
[Fe/H] agree rather well with those of Edvardsson et ah, 
but the overabundance of Na and Al for metal-poor stars 
found in their work is not confirmed. Furthermore, the 
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Galactic evolution of elements not included in Edvards- 
son et al., K, V and Cr, is studied. It is concluded that 
the terms "a elements" and "iron-peak elements" cannot 
be used to indicate production and evolution by specific 
nucleosynthesis processes; each element seems to have a 
unique enrichment history. 

Key words: Nuclear reactions, nucleosynthesis, abun- 
dances - Stars: abundances - Stars: kinematics - Galaxy: 
abundances - Galaxy: evolution - Galaxy: solar neigh- 
bourhood 



1. Introduction 

The chemical abundances of long-lived F and G main se- 
quence stars, combined with kinematical data and ages, 
provide a powerful way to probe the chemical and dy- 
namical evolution of the Galaxy. 

As far as the disk stars are concerned, many general 
trends have been discovered during the past decades. Most 
notable results are correlations of metallicity with age, 
Galactocentric distance, and vertical distance from the 
Galactic plane based on photometric or low-resolution ob- 
servations (e.g. Eggen et al. 1962; Twarog 198C). In addi- 
tion, the abundance patterns for some elements have been 
derived for small samples of stars: oxygen and a elements 
relative to iron vary systematically from overabundances 
at [Fe/H] ~ -1.0 to a solar ratio at [Fe/H] ~ 0.0, while 
most iron-peak elements follow iron for the whole metal- 
licity range of the disk. These results have provided im- 
portant constraints on chemical evolution models for the 
Galactic disk. 

With improved observation and analysis techniques, 
which make it possible to study the Galactic chemical evo- 
lution (GCE) in detail, some old conclusions have, how- 
ever, been challenged and new questions have arisen. Par- 
ticularly important is the detailed abundance analysis of 
189 F and G d warfs w ith -1.1 < [Fe/H] < 0.25 by Ed- 
vardsson et al. ( p.9934 hereafter EAGLNT). The main re- 
sults from this work may be summarized as follows: (1) 
There are no tight relations between age, metallicity and 
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kinematics of disk stars, but substantial dispersions im- 
posed on weak statistical trends. (2) There exists a real 
scatter in the run of [a/Pe] vs. [Fe/H] possibly due to 
the mixture of stars with different origins. The scatter 
seems to increase with decreasing metallicity starting at 
[Fe/H] ~ —0.4. Together with a possible increase in the 
dispersion of TFlsr (the stellar velocity perpendicular to 
the Galactic plane with respect to the Local Standard of 
Rest, LSR) at this point, the result suggests a dual model 
for disk formation. It is, however, unclear if the transition 
at [Fe/H] ~ —0.4 represents the division between the thin 
disk and the thick disk. (3) A group of metal-poor disk 
stars with R m < 7 kpc is found to have larger [a/Fe] values 
than stars with R m > 9 Kpc, indicating a higher star for- 
mation rate (SFR) in the inner disk than that in the outer 
disk. Since essentially all the oldest stars in EAGLNT have 
small i? m , it is, however, difficult to know upon which, 
R m or age, the main dependence of [a/Fe] is. (4) At a 
given age and R m , the scatter in [a/Fe] is negligible while 
[Fe/H] does show a significant scatter. The former implies 
that the products of supernovae of different types are thor- 
oughly mixed into the interstellar medium (ISM) before 
significant star formation occurs. Based on this, the large 
scatter in [Fe/H] may be explained by infall of unpro- 
cessed gas with a characteristic mixing time much longer 
than that of the gas from supernovae of different types. 
(5) The Galactic scatter may be different for individual 
a elements; [Mg/Fe] and [Ti/Fe] show a larger scatter at 
a given metallicity than [Si/Fe] and [Ca/Fe]. It suggests 
that individual a elements may have different origins. (6) 
A new stellar group, rich in Na, Mg, Al, was found among 
the metal-rich disk stars, suggesting additional synthesis 
sources for these elements. 

Given that the study of EAGLNT was based on a lim- 
ited sample of stars with certain selection effects and that 
the analysis technique induced uncertainties in the final 
abundances, some subtle results need further investigation 
before they can provide reliable constraints on theory. For 
example, it is somewhat unclear if the different [a/Fe] at a 
given metallicity between the inner disk and the outer disk 
stars is real and if old disk stars are always located in the 



inner disk. Moreover, recent work by Tomkin et al. (1997) 
argued against the existence of NaMgAl stars. In addition, 
a number of elements, which are highly interesting from 
a nucleosynthetic point of view, were not included in the 
work of EAGLNT. 

The present work, based on a large differently se- 
lected sample of disk stars, aims at exploring and extend- 
ing the results of EAGLNT with improved analysis tech- 
niques. Firstly, we now have more reliable atmospheric 
parameters. The effective temperature is derived from the 
Stromgren b—y color index using a recent infrared- flux cal- 
ibration and the surface gravity is based on the Hipparcos 
parallax. About one hundred iron lines (instead of ~ 30 
in EAGLNT) are used to provide better determinations of 
metallicity and microturbulence. Secondly, the abundance 



calculation is anchored at the most reliable theoretical or 
experimental oscillator strengths presently available in the 
literature. Thirdly, greater numbers of Fen, Sil and Cal 
lines in our study should allow better abundance deter- 
minations, and new elements (K, Sc, V, Cr and Mn) will 
give additional information on Galactic evolution. Lastly, 
the stellar age determination is based on new evolution- 
ary tracks, and the space velocity is derived from more 
reliable distance and proper motion values. 

In the following sections || to ^, we describe the obser- 
vations and methods of analysis in details and present the 
derived abundances, ages and kinematics. The results are 
discussed in Sect. [?] and compared to those of EAGLNT. 
Two elements, Sc and Mn, not included in EAGLNT and 
represented by lines showing significant hyperfine struc- 
ture (HF S) effe cts, are discussed in a separate paper (Nis- 
sen et al. [1999 ) , which includes results for halo stars from 
Nissen & Schuster (|l997j). 



2. Observations 



2.1. Selection of stars 



The stars were selected from the uvby — [3 photometric 
catalogues of Olsen ( 1983| , 1993| ) according to the crite- 



ria of 5800 < T off < 6400 K, 4.0 < \ogg < 4.5 and 
— 1.0 < [Fe/H] < +0.3 with approximately equal numbers 
of stars in every metallicity interval of 0.1 dex. In this se- 
lection, the temperature was determined from the b — y 



index with the calibration of Magain (1987), gravity was 
calculated from the c\ index as described in EAGLNT, 
and metallicity was derived from the mi index using the 
calibrations of Schuster & Nissen (1989). The later rede- 
terminations of the temperature with the calibration of 



Alonso et al. (1996) and the gravity from Hipparcos par- 
allax lead to slight deviations from the selection criteria 
for some stars. 

Based on the above selection, 104 F and G stars were 
observed, but 3 high-rotation (Vsini > 25kms~ 1 ) stars 
and 9 double-line spectroscopic binaries were excluded 
from the sample. Another 11 stars have radial velocity 
dispersions higher than the measurement error of the 
CORAVEL survey. HP 106516 A an d HD 97916 are sus- 
pected binaries (Carney et al. 1994), and HD 25998 and 
HD 206301 are p ossibly variables (e.g. Petit |1990| ; Mor- 
ris & Mutel 1988). These 15 stars (marked in the column 
"Rem" of Table 3 are being carefully used in our study. 
The remaining stars are considered as single stars, but are 
checked for differences in iron abundances between Fel 
and Fe n lines using gravities from Hipparcos parallaxes 
as suggested by Fuhrmann ( |1998| ). As described later, ad- 
ditional 2 stars were excluded during the analysis and thus 
the sample contains 90 stars for the final discussion and 
conclusions. 
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Fig. 1. Examples of spectra obtained with the 2.16m tele- 
scope at Xinglong Station for HD 142373 (T cff = 5920 K, 
\ogg = 4.27, [Fe/H] = -0.39) with a high S/N - 400 and 
HD 106516A (T cff = 6135 K, log 5 = 4.34, [Fe/H] = -0.71) 
with a relatively low S/N ~ 160. 



2.2. Observations and data reduction 

The observations were carried out with the Coude Echelle 
Spectrograph attached to the 2.16m telescope at Beijing 
Astronomical Observatory (Xinglong, P.R. China). The 
detector was a Tek CCD (1024 x 1024 pixels with 24 x 
24 (im 2 each in size). The red arm of the spectrograph with 
a 31.6 grooves/mm grating was used in combination with 
a prism as cross-disperser, providing a good separation 
between the echelle orders. With a 0.5 mm slit (1.1 arcsec), 
the resolving power was of the order of 40 000 in the middle 
focus camera system. 

The program stars were observed during three runs: 
March 21-27, 1997 (56 stars), October 21-23, 1997 (27 
stars) and August 5-13, 1998 (21 stars). The exposure 
time was chosen in order to obtain a signal-to-noise ra- 
tio of at least 150 over the entire spectral range. Most 
bright stars have S/N ~ 200 - 400. Figure shows the 
spectra in the region of the oxygen triplet for two repre- 
sentative stars HD 142373 and HD 106516A. In addition, 
the solar flux spectrum as reflected from the Moon was 
observed with a S/N ~ 250 and used as one of the "stan- 
dard" stars in determining oscillator strengths for some 
lines (see Sect. 4.2). 

The spectra were reduced with standard MIDAS rou- 
tines for order identification, background subtraction, flat- 
field correction, order extraction and wavelength calibra- 
tion. Bias, dark current and scattered light corrections 



Fig. 2. A comparison of equivalent widths measured in 
this work with ESO data in EAGLNT for 5 stars in com- 
mon. The thick line is the linear fit to the points, whereas 
the dashed line is the one-to-one relation. 



are included in the background subtraction. If an early 
B-type star could be observed close to the program stars, 
it was used instead of the flat-field in order to remove in- 
terference fringes more efficiently. The spectrum was then 
normalized by a continuum function determined by fitting 
a spline curve to a set of pre-selected continuum windows 
estimated from the solar atlas. Finally, correction for ra- 
dial velocity shift, measured from at least 20 lines, was 
applied before the measurement of equivalent widths. 

2.3. Equivalent widths and comparison with EAGLNT 

The equivalent widths were measured by three methods: 
direct integration, Gaussian and Voigt function fitting, de- 
pending on which method gave the best fit of the line 
profile. Usually, weak lines are well fitted by a Gaussian, 
whereas stronger lines in which the damping wings con- 
tribute significantly to their equivalent widths need the 
Voigt function to reproduce their profiles. If unblended 
lines are well separated from nearby lines, direct integra- 
tion is the best method. In the case of some intermediate- 
strong lines, weighted averages of Gaussian and Voigt fit- 
ting were adopted. 

The accuracy of the equivalent widths is estimated 
by comparing them to the independent measurements by 
EAGLNT for 25 stars in common. Five of them were 
observed at the ESO Observatory (R ~ 60 000, S/N ~ 
200) and 23 were observed at the McDonald Observatory 
(R ~ 30 000, S/N - 200 - 500). The systematic difference 
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between the two sets of measurements is small and a linear 
least squares fitting gives: 



EWxi = 1.025 (±0.012) SWeso- 
EWxi = 1-083 (±0.006) EW Mc u 



0.89 (±0.56) (mA) 
■0.94 (±0.28) (mA) 



The standard deviations around the two relations are 3.8 
mA (for 129 lines in common with ESO) and 4.3 mA (for 
575 lines in common with McDonald). Given that the er- 
ror of the equivalent widths in EAGLNT is around 2 mA, 
we estimate an rms error of about 3 mA in our equivalent 
widths. As shown by the comparison of our equivalent 
widths with ESO data in Fig. ^, the equivalent widths 
below 50 mA are consistent with the one-to-one relation. 
The deviations for the stronger lines may be due to the 
fact that all lines in EAGLNT were measured by Gaus- 
sian fitting, which leads to an underestimate of equivalent 
widths for intermediate-strong lines because of neglecting 
their wings. We conclude that the Xinglong data may be 
more reliable than the EAGLNT data for lines in the range 
of 50 < EW < 100 mA. 



3. Stellar atmospheric parameters 

3.1. Effective temperature and metallicity 

The effective temperature was determined from the 
Stromgren indices (b — y and ci) and [Fe/H] using the 
calibration of Alonso et al. (1996). If the color excess 
E(b — y), as calculated from the Hp index calibration by 
Olsen ( 1988] ), is larger than 0.01, then a reddening correc- 
tion was applied. 

The metallicity, required in the input for tempera- 
ture and abundance calculation, was first derived from the 
Stromg ren m i index using the calibrations of Schuster & 
Nissen ( 198S ). But the spectroscopic metallicity obtained 
later was used to iterate the whole procedure. 

The errors of the photometric data are a(b—y) = 0.004 
and <t(ci) = 0.008 according to Olsen ( 1993| ). Adopting 
er([Fe/H]) = 0.1 from the spectroscopic analysis, the sta- 
tistical error of T e g is estimated to about ±50 K. Con- 
sidering a possible error of ±50 K in the calibration, the 
error in temperature could reach ±70 K. We do not adopt 
the excitation temperature, determined from a consistent 
abundance derived from Fei lines with different excita- 
tion potentials, because errors induced by incorrect damp- 



ing parameters (Ryan 1998) or non-LTE effects can be 
strongly dependent on excitation potential, leading to an 
error in effective temperature as high as 100 K. 

3.2. Gravity 

In most works, gravities are determined from the abun- 
dance analysis by requiring that Fei and Fell lines give 
the same iron abundance. But it is well known that the 
derivation of iron abundance from Fe I and Fe II lines may 



be affected by many factors such as unreliable oscilla- 
tor strengths, possible non-LTE effects and uncertainties 
in the temperature structure of the model atmospheres. 
From the Hipparcos parallaxes, we can determine more 
reliable gravities using the relations: 
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M bol = V^ + BC + SlogTr + S, 



(2) 



where, M. is the stellar mass, Mbol the absolute bolometric 
magnitude, V the visual magnitude, BC the bolometric 
correction, and 7r the parallax. 

The parallax is taken from the Hipparcos Satellite ob- 
servations (ESA 1997 ). For most program stars, the rela- 
tive error in the parallax is of the order of 5%. Only two 
stars in our sample have errors larger than 10%. From 
these accurate parallaxes, stellar distances and absolute 
magnitudes were obtained. Note, however, that our sam- 
ple includes some binaries, for which the absolute magni- 
tude from the Hipparcos parallax could be significantly in 
error. An offset of —0.75 mag will be introduced for a bi- 
nary with equal components through the visual magnitude 
in Eq. (0). Thus, we also calculated absolute magnitudes 
from the photometric indices f3 and c\ using the relations 
found by EAGLNT. Although the absolute magnitude of 
a binary derived by the photometric method is also not 
very accurate due to different spectral types and thus dif- 
ferent flux distributions of the components, it may be bet- 
ter than the value from the parallax method. Hence, for 
a few stars with large differences in absolute magnitudes 
between the photometry and parallax determination, we 
adopt the photometric values. 

The bolometric correctio n was interpolated from the 
new BC grids of Alonso et al. ( 1995| ) determined from line- 
blanketed flux distributions of ATLAS9 models. It is noted 
that the zero-point of the bolometric correction adopted 
by Alonso et al., BCq = —0.12, is not consistent with 
the bolometric magnitude of the Sun, Af D oi,o— 4.75, re- 
cently recommended by the IAU ( |1999| ). But the gravity 
determination from the Eq. (|l|) only depends on the M DO i 
difference between the stars and the Sun and thus the 
zero-point is irrelevant. 

The derivation of mass is described in Sect. |[ The 
estimated error of 0.05 M Q in mass corresponds to an 
error of 0.03 dex in gravity, while errors of 0.05 mag in 
BC and 70 K in temperature each leads to an uncertainty 
of 0.02 dex in logg. The largest uncertainty of the gravity 
comes from the parallax. A typical relative error of 5% 
corresponds to an error of 0.04 dex in logg. In total, the 
error of logg is less than 0.10 dex. 

The surface gravity was also estimated from the 
Balmer discontinuity index c\ as described in EAGLNT. 
We find a small systematical shift (about 0.1 dex) be- 
tween the two sets of logg, with lower gravities from 
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the parallaxes. There is no corresponding shift between 
My(par) and My (phot). The mean deviation is 0.03 mag 
only, which indicates that the systematic deviation in log<? 
comes from the gravity calibration in EAGLNT. 

3.3. Microturbulence 

The microturbulence, £t, was determined from the abun- 
dance analysis by requiring a zero slope of [Fe/H] vs. EW. 
The large number of Fei lines in this study enables us 
to choose a set of lines with accurate oscillator strengths, 
similar excitation potentials (xiow > 4.0 eV) and a large 
range of equivalent widths (10 - 100 mA) for the determi- 
nation. With this selection, we hope to reduce the errors 
from oscillator strengths and potential non-LTE effects for 
Fe I lines with low excitation potentials. The error of the 
microturbulence is about 0.3 kms . 

The relation of £t as a function of T e g and logg de- 
rived by EAGLNT corresponds to about 0.3 kms -1 lower 
values than those derived from our spectroscopic analy- 
sis. No obvious dependence of the difference on tempera- 
ture, gravity and metallicity can be found. In particular, 
the value for the Sun in our work is 1.44 kms -1 , also 
0.3 kms -1 higher than the value of 1.15 found from the 
EAGLNT relation. The difference in & between EAGLNT 
and the present work is probably related to the difference 
in equivalent widths of intermediate-strong lines discussed 
in Sect. 2.3, EAGLNT measured these lines by fitting a 
Gaussian function and hence underestimated their equiv- 
alent widths, leading to a lower microturbulence. 

Finally, given that the atmospheric parameters were 
not determined independently, the whole procedure of de- 
riving T ff, logg, [Fe/H] and £ t was iterated to consis- 
tency. The atmospheric parameters of 90 stars are pre- 
sented in Table 3. The uncertainties of the parameters 
are: er(T cff ) = 70 K, aQogg) = 0.1, cr([Fe/H]) = 0.1, and 
cr(&) = 0.3 kms -1 . 



4. Atomic line data 

4-1. Spectral lines 

All unblended lines with symmetric profiles having equiv- 
alent widths larger than 20 mA in the solar atlas (Moore 
et al. 1966 ) were cautiously selected. The equivalent width 
limit ensures that lines are not disappearing in the most 
metal-poor disk stars at [Fe/H] = —1.0. Given that very 
weak lines would lead to an increase of random errors in 
the abundance determination and that too strong lines 
are very sensitive to damping constants, only weak and 
intermediate-strong lines with 3 < EW < 100 mA in the 
stellar spectra were adopted in our abundance analysis 
except for potassium, for which only one line (K I A7699) 
with an equivalent width range of 50 - 190 mA, is avail- 
able. 



4-2. Oscillator strengths 

Due to the large number of measurable lines in the spec- 
tra, Fe I lines were used for microturbulence determina- 
tion and temperature consistency check. Hence, careful 
selection of oscillator strengths for them is of particular 
importance. Of many experimental or theoretical calcula- 
tions of oscillator strengths for Fe I lines, only three sources 
with precise gf values were chosen. They are: Blackwell 
et al. ( h982b| ; |l982d) , O 'Bria n et a l. ( |1991| ) and Bard & 
Kock ( |199lD or Bard et al. (|l994j) . The agreements be- 
tween these sources are very satisfactory, and thus mean 
logy/ values were adopted if oscillator strengths are avail- 
able in more than one of the three sources. A few oscilla- 
tor strengths with large differences between these sources 
were excluded. 

References for other elements are: Oi (Lambert 1978) , 
Nai (Lambert & W arner |196S| ), Mg i (Ch ang |1990[ ), Al l 



(Lambert fe Wa rner |i~968| ) , Sil (Garz |l973j), K i ( Lambert 
& Warner [1968D Cai (Sm ith & Raggett |1 98 1| ), Ti I (B lack- 
well et al. |!982a| |l986a| ), Vi (Whaling et al. |1985[ ), Cri 
(Blackwell et al. 1986b| ), Fen ( Biemo nt et al. 1991; Han- 
nafor d et al. |1992| ), Nil (Kostyk |i982] Wickliffe & Lawler 
|1997[ ) and Ban (Wiese & Martin |198C| ). 

These experimental or theoretical gf values were in- 
spected to see if they give reliable abundances by evalu- 
ating the deviation between the abundance derived from 
a given line and the mean abundance from all lines of the 
same species. A significant mean deviation in the same 
direction for all stars (excluding the suspected binaries) 
were used to correct the gf value. Lines with large devi- 
ations in different directions for different stars were dis- 
carded from the experimental or theoretical logg/ list. The 
two sets of values are presented in columns "abs." and 
"cor." in Table 4. Based on these gf values, we derived 
abundances for 10 "standard" stars: HD 60319, HD 49732, 
HD 215257, HD 58551, HD 101676, HD 76349, HD 58855, 
HD 22484, the Sun and HD 34411 A (in metallicity order). 
Oscillator strengths for lines with unknown gf values were 
then determined from an inverse abundance analysis of 
the above 10 stars, which are distributed in the metal- 
licity range —0.9 < [Fe/H] < 0.0 dex and were observed 
at high S/N — 250 - 400. Generally, the gf values for a 
given line from different "standard" stars agree well, and 
a mean value (given in the column "dif" in Table 4 was 
thus adopted. 

4-3. Empirical enhancement factor 

It has been recognized for a long time that line broading 



derived from Unsold's (1955) approximation to Van der 
Waals interaction is too weak, and an enhancement factor, 
E~ f , should be applied to the damping parameter, 7. Usu- 
ally, enhancement factors of Fe 1 lines with excitation po- 
tential of the lower energy level (xiow) less than 2.6 eV are 
taken from the empirical calibration by Simmons & Black- 
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well fll982| ). For Fei lines with xiow > 2.6 eV, E 1 = 1.4 
is generally used in abundance analysis. Recently, Anstee 
&; O'Mara's ( 1995| ) computed the broadening cross sec- 



5.2.1. Errors from equivalent widths and atomic data 



The comparison of equivalent widths in Sect. 2.3 indicates 



tions for s-p and p-s transitions and found that Ey should 
be 2.5 ~ 3.2 for lines with xiow > 3.0 cV, whereas lines 
with xiow < 2.6 eV have broadening cross sections more 
consistent with Simmons & Blackwcll's (1982) work. 

Following EAGLNT and other works, we adopted Sim- 
mons & Blackwell's ( |1982| ) E 1 for Fei lines with xiow < 2.6 
eV and E~ — 1.4 for the remaining Fe I lines, while a value 
of 2.5 was applied for Fen lines as suggested by Holweger 



et al. (1990). Enhancement factors for Nal, Si I, Cai and 



Ban were taken from EAGLNT (see references therein). 
Finally, a value of 1.5 was adopted for the Ki, Til and 
Vi lines considering their low excitation potentials, and 
a factor of 2.5 was applied to the remaining elements fol- 
lowing Mackle et al. ( 1975 ). The effects of changing these 
value s by 50% on the derived abundances are discussed in 
Sect. ~~ 



5.2 



The atomic line data are given in Table 4. 



5. Abundances and their uncertainties 

5.1. Model atmospheres and abundance calculations 

The abundance analysis is based on a grid of flux con- 
stant, homogeneous, LTE model atmospheres, kindly sup- 
plied by Bengt Edvardsson (Uppsala). The models were 
computed with the MARCS code using the updated con- 



tinuous opacities by Asplund et al. (1997) including UV 
line blanketing by millions of absorption lines and many 
molecular lines. 

The abundance was calculated with the program 
EQWIDTH (also made available from the stellar atmo- 
spheric group in Uppsala) by requiring that the calcu- 
lated equivalent width from the model should match the 
observed value. The calculation includes natural broad- 
ening, thermal broadening, van der Waals damping, and 
the microturbulent Doppler broadening. The mean abun- 
dance was derived from all available lines by giving equal 
weight to each line. Finally, solar abundances, calculated 
from the Moon spectrum, were used to derive stellar abun- 
dances relative to solar values (Table 5). Such differen- 
tial abundances are generally more reliably than absolute 
abundances because many systematic errors nearly cancel 
out. 

5.2. Uncertainties of abundances 

There are two kinds of uncertainties in the abundance 
determination: one acts on individual lines, and includes 
random errors of equivalent widths, oscillator strengths, 
and damping constants; another acts on the whole set of 
lines with the main uncertainties coming from the atmo- 
spheric parameters. 



that the typical uncertainty of the equivalent width is 
about 3 mA, which leads to an error of about 0.06 dex 
in the elemental ratio X/H derived from a single line with 
an equivalent width around 50 mA. For an element rep- 
resented by N lines, the error is decreased by a factor 
y/W. In this way, the errors from equivalent widths were 
estimated for elements with only one or a few lines. Al- 
ternatively, the scatter of the deduced abundances from 
a large number of lines with reliable oscillator strengths 
gives another estimate of the uncertainty from equivalent 
widths. With over 100 Fei lines for most stars, the scatter 
varies somewhat from star to star with a mean value of 
0.07 dex, corresponding to an error of 0.007 dex in [Fe/H]. 
Other elements with significant numbers of lines, such as 
Ca, Ni and Si, have even smaller mean line-to-line scatters. 

The uncertainties in atomic data are more difficult 
to evaluate. But any error in the differential abundance 
caused by errors in the gf values is nearly excluded due 
to the correction of some experimental or theoretical gf 
values and the adoption of mean gf values from 10 "stan- 
dard" stars. Concerning the uncertainties in the damping 
constants, we have estimated their effects by increasing 
the adopted enhancement factors by 50%. The microtur- 
bulence was accordingly adjusted because of the coupling 
between the two parameters. The net effect on the differ- 
ential abundances with respect to the Sun is rather small 
as seen from Table [l| 

5.2.2. Consistency check of atmospheric parameters 

As a check of the photometric temperature, the derived 
iron abundance from individual Fe I lines was studied as a 
function of the excitation potential. To reduce the influ- 
ence of microturbulence, only lines with equivalent widths 
less than 70 mA were included. A linear least squares 
fit to the abundance derived from each line vs. low ex- 
citation potential determines the slope in the relation 
[Fe/H] = a+b-xiow- The mean slope coefficient for all stars 
is b = 0.004 ± 0.013. There is only a very small (if any) 
dependence of b on effective temperature, surface gravity 
or metallicity. A suspected binary, HD 15814, has a very 
deviating slope coefficient (b = —0.056) and is excluded 
from further analysis. 

The agreement of iron abundances derived from Fei 
and Fe n lines is satisfactory when gravities based on Hip- 
parcos parallaxes are used (see Fig. [|). The deviation is 
less than 0.1 dex for most stars with a mean value of 
-0.009 ±0.07 dex. From [Fe/H] = 0.0 to [Fe/H] = -0.5, 
the mean deviation ([Fe/H] n — [Fe/H]j) seems, however, 
to increase by about 0.1 dex in rough agreement wit h pre- 
dictions from non-LTE computations (see Sect. 5.3). 

The deviation in iron abundances based on Fe I and 
Fe II abundance provides a way to identify binaries and 
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Fig. 3. Difference in iron abundances derived from Fel 
and Fen lines vs. [Fe/H] with suspected binaries marked 
by a square around the filled circles. 



to estimate the influence of the component on the pri- 
mary. The suspected binaries are marked with an addi- 
tional square around the filled circles in Fig. 0. It shows 
that there is no significant influence from the component 
for these binaries except in the case of HD 15814, which 
was already excluded on the basis of it's ^-coefficient in the 
excitation equilibrium of Fel lines. Thus, the other possi- 
ble binaries were included in our analysis. It is, however, 
surprising that HD 186257 show a higher iron abundance 
based on Fe n lines than that from Fe I lines with a devia- 
tion as large as 0.28 dex. We discard this star in the final 
analysis and thus have 90 stars left in our sample. 

5.2.3. Errors in resulting abundances 

Table [l] shows the effects on the derived abundances of a 
change by 70 K in effective temperature, 0.1 dex in gravity, 
0.1 dex in metallicity, and 0.3 kms -1 in microturbulence, 
along with errors from equivalent widths and enhancement 
factors, for two representative stars. 

It is seen that the relative abundances with respect to 
iron are quite insensitive to variations of the atmospheric 
parameters. One exception is [O/Fe] due to the well known 
fact that the oxygen abundance derived from the infrared 
triplet has an opposite dependence on temperature to that 
of the iron abundance. After rescaling of our oxygen abun- 
dances to results from the forbidden line at A6300 (see 
next section) the error is somewhat reduced. Therefore, 
the error for [O /Fe] in Table |l] might be overestimated. 

In all, the uncertainties of the atmospheric parameters 
give errors of less than 0.06 dex in the resulting [Fe/H] 
values and less than 0.04 dex in the relative abundance 
ratios. For an elemental abundance derived from many 
lines, this is the dominant error, while for an abundance 
derived from a few lines, the uncertainty in the equivalent 



Table 1. Abundance errors. The last column gives the 
total error assuming that the individual errors are uncor- 
rected 
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widths may be more significant. Note that the uncertain- 
ties of equivalent widths for V and Cr (possibly also Ti) 
might be underestimated given that their lines are gener- 
ally weak in this work. In addition, with only one strong 
line for the K abundance determination, the errors from 
equivalent widths, microturbulence and atomic line data 
are relatively large. 

Lastly, we have explored the HFS effect on one All 
line at A6698 , one Mgi line at A5711, and two Ball lines 
at A6141 and A6496. The HFS data are taken from three 
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sources: Bichl ( |1976D for Al, Steffen fll985| ) for Mg and 
Frangois ( 1996j ) for Ba. The results indicate that the HFS 
effects are very small for all these lines with a value less 
than 0.01 dex. 



5.3. Non-LTE effects and inhomogeneous models 

The assumption of LTE and the use of homogeneous 
model atmospheres may introduce systematic errors, es- 
pecially on the slope of various abundance ratios [X/Fe] 
vs. [Fe/H]. These problems were discussed at quite some 
length by EAGLNT. Here we add some remarks based on 
recent non-LTE studies and computations of 3D hydrody- 
namical model atmospheres. 

Based on a number of studies, EAGLNT concluded 
that the maximum non-LTE correction of [Fe/H], as de- 
rived from Fel lines, is 0.05 to 0.1 dex for metal-poor F 



and G disk dwarfs. Recently, Thevenin & Idiart ( pL999[ ) 
computed non-LTE corrections on the order of 0.1 to 0.2 
dex at [Fe/H] = —1.0. Fig. || suggests that the maxim- 
ium correction to [Fe/H] derived from Fel lines is around 
0.1 dex, but we emphasize that this empirical check may 
depend on the adopted T c g calibration as a function of 
[Fe/H] . 

The oxygen infrared triplet lines are suspected to be 
affected by non-LTE formation, because they give system- 
atically higher abundances than forbidden lines. Recent 
work by Reetz (1999) indicates that non-LTE effects are 
insignificant (< 0.05 dex) for metal-poor and cool stars, 
but become important for warm and metal-rich stars. For 
stars with [Fe/H] > —0.5 and T c g > 6000 K in our sample, 
non-LTE effects could reduce the oxygen abundances by 
0.1-0.2 dex. For this reason, we use Eq. (11) of EAGLNT 
to scale the oxygen abundances derived from the infrared 
triplet to those derived by Nissen & Edvardsson (1992) 
from the forbidden [O i] A6300. 

The two weak Nai lines (A6154 and A6160) used for 
our Na abundance determinations, are only marginally af- 
fected by deviations from LTE formation (Baumiillcr ct al. 
1998) ) . The situation for Al may, however, be different. The 



non-LTE analysis by Baumiiller & Gehren ( |1997 ) of one 
of the All lines used in the present work (A6698) leads to 
about 0.15 dex higher Al abundances for the metal-poor 
disk dwarfs than those calculated from LTE. No non-LTE 
study for the other two lines used in the present work 
is available. We find, however, that the derived Al abun- 
dances depend on T e g with lower [Al/Fe] for higher tem- 
perature stars. This may be due to the neglect of non-LTE 
effects in our work. Hence, we suspect that the trend of 
[Al/Fe] vs. [Fe/H] could be seriously affected by non-LTE 
effects. 

The recent non-LTE analysis of neutral magnesium in 
the solar atmosphere by Zhao et al. ( 1998j ) and in metal- 
poor stars by Zhao & Gehren (1999) leads to non-LTE 



corrections of 0.05 dex for the Sun and 0.10 dex for a 
[Fe/H] — — 1.0 dwarf, when the abundance of Mg is de- 



rived from the A5711 Mgi line. Similar corrections are ob- 
tained for some of the other lines used in the present work. 
Hence, we conclude that the derived trend of [Mg/Fe] vs. 
[Fe/H] is not significantly affected by non-LTE. 

The line-profile analysis of the K I resonance line at 
A7699 by Takeda et al. ( ]1996| ) shows that the non-LTE 
correction is —0.4 dex for the Sun and —0.7 dex for Pro- 
cyon. There are no computations for metal-poor stars, but 
given the very large corrections for the Sun and Procyon 
one may expect that the slope of [K/Fe] vs. [Fe/H] could 
be seriously affected by differential non-LTE effects be- 
tween the Sun and metal-poor stars. 

The non-LTE study of Ba lines by Mashonkina et al. 
( 1999| ), which includes two of our three Ban lines (A5853 
and A6496), give rather small corrections (< 0.10 dex) to 
the LTE abundances, and the corrections are very similar 
for solar metallicity and [Fe/H] ~ —1.0 dwarfs. Hence, 
[Ba/Fe] is not affected significantly. 

In addition to possible non-LTE effects, the derived 
abundances may also be affected by the representation 
of the stellar atmospheres by plane-parallel, homogeneous 
models. The recent 3D hydrodynamical model atmo- 
spheres of metal-poor stars by Asplund et al. ( |1999[ ) have 
substantial lower temperatures in the upper photosphere 
than ID models due to the dominance of adiabatic cool- 
ing over radiative heating. Consequently, the iron abun- 
dance derived from Fel lines in a star like HD 84937 
(T cff ~ 6300 K, \ogg ~ 4.0 and [Fe/H] ~ -2.3) is 0.4 dex 
lower than the value based on a ID model. Although the 
effect will be smaller in a [Fe/H] ~ —1.0 star, and the 
derived abundance ratios are not so sensitive to the tem- 
perature structure of the model, we clearly have to worry 
about this problem. 



5.4- Abundance comparison of this work with EAGLNT 

A comparison in abundances between this work and 
EAGLNT for the 25 stars in common provides an inde- 
pendent estimate of the errors of the derived abundances. 
The results are summarized in Table ^. 

The agreement in iron abundance derived from Fel 
lines is satisfactory with deviations within ±0.1 dex for 
the 25 common stars. These small deviations are mainly 
explained by different temperatures given the fact that 
the abundance differences increase with temperature de- 
viations between the two works. The rms deviation in iron 
abundance derived from Fe II lines are slightly larger than 
that from Fe I lines. The usage of different gravities partly 
explain this. But the small line-to-line scatter from 8 Fe n 
lines in our work indicates a more reliable abundance than 
that of EAGLNT who used 2 Fen lines only. 

Our oxygen abundances are systematically higher by 
0.15 dex than those of EAGLNT for 5 common stars. 
Clearly, the temperature deviation is the main reason. The 
systematically lower value of 70 K in our work increases 
[O/Fe] by 0.10 dex (see Table @). 
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Table 2. Mean abundance differences (this 
work— EAGLNT) and standard deviations. N is the 
number of stars, for which a comparison was possible 
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[Al/Fe] 
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[Si/Fe] 


-0.012 
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23 




[Ca/Fe] 


0.055 


0.039 


25 




[Ti/Fe] 


-0.093 


0.099 


24 




[Ni/Fe] 


-0.008 


0.045 


25 


6 - 


[Ba/Fe] 


0.068 


0.081 


25 





The mean abundance differences for Mg, Al, Si, Ca 
and Ni between the two works are hardly significant. The 
systematical differences (this work - EAGLNT) of —0.08 
dex for [Na/Fc] and [Ti/Fe] and +0.07 dcx for [Ba/Fe] 
are difficult to explain, but we note that when the abun- 
dances are based on a few lines only, a systematic offset 
of the stars relative to the Sun may occur simply because 
of errors in the solar equivalent widths. 

6. Stellar masses, ages and kinematics 

6.1. Masses and ages 



As described in Sect. 3.2, the stellar mass is required in the 
determination of the gravity from the Hipparcos parallax. 
With the derived temperature and absolute magnitude, 
the mass was estimated from the stellar position in the 
My — logT c ff diagram (see Fig. ||) by interp olatin g in the 
evolutionary tracks of VandenBerg et al. ( 1999| ), which 
are distributed in metallicity with a step of ~ 0.1 dex. 
These new tracks are based on the recent OPAL opacities 
(Rogers & Iglesias 1992) using a varying helium abun- 
dance with [a/Fe] = 0.30 for [Fe/H] < -0.3 and a con- 
stant helium abundance (Y = 0.2715) without a element 
enhancement for [Fe/H] > —0.2. Fig. ^ shows the position 
of our program stars with —0.77 < [Fe/H] < —0.66 com- 
pared to the evolutionary tracks of Z = 0.004 ([Fe/H] = 
—0.71). The errors in T e ff, My, and [Fe/H] translate to an 
error of 0.06-M© in the mass. 

Stellar age is an important parameter when study- 
ing the chemical evolution of the Galaxy as a function 
of time. Specifically, the age is useful in order to inter- 
pret abundance ratios as a function of metallicity. In this 
work, the stellar age was obtained simultaneously with the 
mass from interpolation in the evolutionary tracks of Van- 
denBerg et al. ( pJ99| ). It was checked that practically the 
same age is derived from the corresponding isochrones. As 
an example, a set of stars are compared to isochrones in 
Fig. |j| The error of the age due to the uncertainties of 




Fig. 4. The positions of our program stars with —0.77 < 
[Fe/H] < —0.66 co mpar ed to the evolutionary tracks of 
VandenBerg et al. (|1999[) with [Fe/H] = -0.71. 



T cff , M v , and [Fe/H] is about 15% (cr(logr) = 0.07) ex- 
cept for a few stars, which have relatively large errors of 
the Hipparcos parallaxes. 



3.5 - 



4.5 




5.5 - 



6.0 



Fig. 5. The positions of our program stars with —0.77 < 
[Fe/H] < —0.66 compared to the isochrones of Vanden- 
Berg et al. ( |1999D with [Fe/H] = -0.71. 
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6.2. Kinematics 

Stars presently near the Sun may come from a wide range 
of Galactic locations. Information on their origin will help 
us to understand their abundance ratios. Therefore, stel- 
lar space velocity, as a clue to the origin of a star in the 
Galaxy, is very interesting. 

The accurate distance and proper motion available 
in the Hipparcos Catalogue (ESA 1997), combined with 
stellar radial velocity, make it possible to derive a reli- 
able space velocity. Radial velocities from the CORAVEL 
survey for 53 stars were kindly made available by Nord- 
strom (Copenhagen) before publication. These velocities 
are compared with our values derived from the Doppler 
shift of spectral lines. A linear least squares fit for 40 stars 
(excluding the suspected binaries) gives: 

RV = 0.997 (±0.002) AVcoravel + 0.26 (±0.12) kms -1 

The rms scatter around the relation is 0.72 kms -1 , show- 
ing that our radial velocities are as accurate as 0.5 kms -1 . 
Hence, our values are adopted for stars not included in the 
CORAVEL survey. 

The calculation of the space velocity with respect to 
the Sun is based on the method presented by Johnson 



& Soderblom (1987). The correction of space velocity to 
the Local Standard of Rest is based on a solar motion, 
(-10.0,+5.2,+7.2) kms -1 in (U,V,W)[ L as d erived from 
Hipparcos data by Dehnen & Binney ( |1998| ). The error 
in the space velocity arising from the uncertainties of dis- 
tance, proper motion and radial velocity is very small with 
a value of about ±1 kms -1 . 

The ages and space velocities derived in the present 
work are generally consistent with EAGLNT. But the 
more accurate absolute magnitude, as well as the new 
set of theoretical isochrones, in our study should give 
more reliable ages than those determined by EAGLNT 
based on the photometric absolut e mag nitude and the old 
isochrones of VandenBerg & Bell (1985). This situation is 
also true for space velocities with our results based on dis- 
tances and proper motions now available from Hipparcos. 

7. Results and discussion 

7.1. Relations between abundances, kinematics and ages 

The observed trends between abundances, kinematics and 
ages are the most important information for theories of 
Galactic evolution. Especially EAGLNT have provided 
many new results on this issue. For example, the substan- 
tial dispersion in the AMR found by EAGLNT argues 
against the assumption of chemical homogeneity adopted 
in many chemical evolution models. It is, however, impor- 
tant to test the results of EAGLNT for a different sample 
of disk stars. Based on more reliable ages and kinematics, 
the present study makes such an investigation. 

1 In the present work, U is defined to be positive in the an- 
ticentric direction. 
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Fig. 6. The abundances of a elements, Fe and Ba as a 
function of logarithmic age. The symbols are shown in 



7.1.1. Age-metallicity relation in the disk 

Fig. |6| shows the age-metallicity relations for a, iron^] and 
barium elements, where a represents the mean abundance 
of Mg, Si, Ca and Ti. Generally, there is a loose correlation 
between age and abundance. Stars younger than 5 Gyr 
(logrg < 0.7) are more metal-rich than [Fe/H] ~ —0.3, 
and stars with [Fe/H] < —0.5 are not younger than 6- 
7 Gyr (logrg > 0.8). The deviating young halo star 
HD 97916 (indicated by an asterisk in Fig. ||) is discussed 
in Sect. |. 

The correlation between age and abundance is, how- 
ever, seriously distorted by a considerable scatter. Stars 
with solar metallicity have an age spread as large as 10 
Gyr, and coeval stars at 10 Gyr show metallicity differ- 
ences as high as 0.8 dex. Such a dispersion cannot be 
explained by either the abundance error (< 0.1 dex) or 
the age uncertainty (~ 15%) in the AMR. This is an im- 



2 Here and in the following sections and figures, the iron 
abundance is the mean abundance derived from all Fel and 
Fell lines with equal weight to each line. 
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portant constraint on GCE models, which must reproduce 
both the weak correlation and the substantial dispersion. 

It is seen from Fig. || that Ba has the steepest slope in 
the AMR, Fe has intermediate slope, and the a elements 
show only a very weak trend with [Fe/H]. This was also 
found by EAGLNT and is consistent with nucleosynthesis 
theory that suggests that the main synthesis sites of Ba, Fe 
and a elements are AGB stars (1-3 M Q ), SNe la (6-8 M Q ) 
and SNe II (> 8M©), respectively. Due to their longer 
lifetime, lower mass stars contribute to the enrichment of 
the Galaxy at a later epoch, i.e. after massive stars have 
been polluting their products into the ISM. Hence, the Ba 
abundance is relatively low in the beginning of the disk 
evolution and increases quickly in the late stage, leading 
to a steeper slope. 

It is interesting that there is a hint of a smaller metal- 
licity spread for young stars with logrg = 0.4 — 0.8 in 
this work than in EAGLNT, while the spread is similar 
for old stars. If we are not misled by our sample (less 
young stars than in the EAGLNT sample and a lack of 
stars with logr 9 < 0.4), it seems that there are metal- rich 
stars at any time in the solar neighbourhood while metal- 
poor stars are always old. Another interesting feature for 
the young stars is that [Ba/H] has a smaller metallicity 
spread than [Fe/H] and [a/H]. This could be due to the 
dependence of elemental yield on the progenitor's mass. 
Ba is produced by AGB stars with a rather small mass 
range of 1-3 Mq , while Fe and a elements are synthesized 
by SNe having a mass range ~6-30 Mq. 

7.1.2. Stellar kinematics as functions of age and 
metallicity 

The study of the dispersion in kinematical parameters as 
a function of Galactic time is more interesting than the 
kinematical data alone, because any abrupt increase in 
dispersion may indicate special Galactic processes occur- 
ring during the evolution. Generally, dispersions in Vlsr, 
IVlsr and total velocity increase with stellar age. We have 
not enough stars at ~ 2.5 Gyr and 10 Gyr to confirm 
the abrupt increases in the Wlsr dispersion found by 
EAGLNT at these ages. Instead, our data seems to in- 
dicate that the kinematical dispersion (possibly also the 
metallicity) is fairly constant for stars younger than 5 
Gyr (logrg = 0.7), but it increases with age for stars 
with logrg > 0.7. Coincidentally, 5 Gyr corresponds to 
[Fe/H] ~ -0.4 dex, the metallicity where EAGLNT sug- 
gested an abundance transition related to a dual forma- 
tion of the Galactic disk. The abundance transition at 
[Fe/H] ~ —0.4 dex is confirmed by our data, but the in- 
crease of the Wlsr dispersion at [Fe/H] ~ —0.4 found by 
EAGLNT is less obvious in our data. 

When the velocity component in the direction of 
Galactic rotation, Vlsr, is investigated as a function of 
the metallicity (see Fig. 0), we find that there are two 
subpopulations for [Fe/H] < —0.6 with positive Vlsr in 
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Fig. 7. Vlsr vs. [Fe/H] with different symbols showing 
three groups of stars. The asterisk indicates the halo star 
HD 97916 discussed in Sect. H. 



group A and negative Vlsr in group C, while stars with 
[Fe/H] > -0.6 have V LS r around Vlsr = -lOkms" 1 
(group B). The pattern persists when other elements 
are substituted for Fe. As shown in Edvardsson et al. 
( |l993h| ), there is a tight correlation between Vlsr and 
the mean Galactocentric distance in the stellar orbit, R m . 
Hence, we can trace the metallicity at different Galacto- 
centric distances assuming that R m is a reasonable es- 
timator of the radius of the star's original orbit. Note, 
however, that the lower metallicity toward the Galactic 
center (Vlsr ^ — 50kms _1 ) for group C stars may be 
due to their large ages. Excluding these stars, a trend of 
decreasing metallicity with increasing Vlsr for stars with 
similar age is found, which indicates a radial abundance 
gradient in the disk, and thus suggests a faster evolution 
in the inner disk than the outer. This is compatible with 
a higher SFR, due to the higher density, in the inner disk. 

There are two possibilities to explain the stars in group 
C. One is anchored to the fact that the oldest stars (> 10 
Gyr) in our sample have the lowest Vlsr, i-e. the smallest 
R m , indicating that the Galaxy did not extend to the Sun 
at 10 Gyr ago according to an inside-out formation process 
of the Galaxy. The other is that these stars come from the 
thick disk, which is older and more metal-poor than the 
thin disk. 

7.2. Relative abundances 

The general trends of elemental abundance with respect 
to iron as a function of metallicity, age and kinematics 
are to be studied in connection with Galactic evolution 
models and nucleosynthesis theory. The main results are 
shown in Fig. || and will be discussed together with those 
of EAGLNT. 
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7.2.1. Oxygen and magnesium 



7.2.2. Silicon, calcium and titanium 



In agreement with most works, [O/Fe] shows a tendency 
to decrease constantly with increasing metallicity for disk 
stars. As oxygen is only produced in the massive progeni- 
tors of SNe II, lb and Ic, it is mainly build up at early times 
of the Galaxy, leading to an overabundance of oxygen in 
halo stars. The [O/Fe] ratio gradually decline in the disk 
stars when iron is produced by the long-lived SNe la. The 
time delay of SNe la relative to SNe II is responsible for 
the continuous decrease of oxygen in disk stars. The ten- 
dency for [O/Fe] to continue to decrease at [Fe/H] > —0.3 
argues for an increasing ratio of SNe la to SNe II also at 
the later stages of the disk evolution. 

In general, the relation of [O/Fe] vs. Vlsr reflects the 
variation of [Fe/H] with Vlsr (see Fig. 0) . [O /Fe] decreases 
with increasing Vlsr for stars with Vlsr < and slowly 
increase with further larger Vlsr- Considering their simi- 
lar ages, the decreasing [O/Fe] from group A to group B 
stars may be attributed to the increasing Vlsr, whereas 
the higher [O/Fe] of group C is due to an older age. 

The magnesium abundance shows a decreasing trend 
with increasing metallicity like oxygen for [Fe/H] < —0.3 
but it tends to flatten out for higher metalicities. Given 
that magnesium is theoretically predicted to be formed 
only in SNe II, the similar decreasing trend as oxygen 
is easily understood, but the flat [Mg/Fe] towards higher 
metallicities than [Fe/H] > —0.3 is unexpected. It seems 
that SNe II are not the only source for Mg. Perhaps SNe 
la also contribute to the enrichment of Mg during disk 
evolution. 

The flat trend of [Mg/Fe] vs. [Fe/H] for [Fe/H] > -0.3 
is also evident from the data of EAGLNT if the high 
Mg/Fe ratios of their NaMgAl stars are re duced to a so- 
lar ratio as found by Tomkin et al. (1997). Feltzing and 
Gustafsson ( 1999D also find [Mg/Fe] to be independent of 
metallicity for their more metal-rich stars although the 
scatter is large. 

With more magnesium lines in the present study, we 
get a similar scatter of [Mg/Fe] as EAGLNT. The scatter 
is slightly larger than that of oxygen and in particular 
much larger than those of Si and Ca. Although we do 
not find a large line-to-line scatter in the Mg abundance 
determination, it is still unclear if the scatter in [Mg/Fe] 
is cosmic. Only 3 Mgi lines are available for most stars 
while Si and Ca are represented by 20-30 lines. There is 
no obvious evidence showing the scatter to be an effect of 
different Vlsr- Nor do we find a clear separation of thick 
disk stars from thin disk stars in the diagram o f [Mg /Fe] 
vs. [Mg/H], as has been found by Fuhrmann ( |l998|) . It 
seems that neither observation nor theory is satisfactory 
for Mg. 



Like magnesium, [Si/Fe] and [Ca/Fe] decrease with in- 
creasing metallicity for [Fe/H] < —0.4 and then flatten 
out with further increasing [Fe/H]. The result is in agree- 
ment with EAGLNT, who found a "kink" at [Fe/H] = 
—0.3 ~ —0.2. But [Ca/Fe] possibly continues to decrease 
for [Fe/H] > —0.4 based on our data. The suspicion 
that Si is about 0.05 dex overabundant relative to Ca for 
[Fe/H] > —0.2 and the possible upturn of silicon at higher 
metallicity in EAGLNT are not supported by our data. 

Both Si and Ca have a very small star-to-star scatter 
(0.03 dex) at a given metallicity for thin disk stars. The 
scatter is slightly larger among the thick disk stars. Since 
the scatter corresponds to the expected error from the 
analysis, we conclude that the Galactic scatter for [Si/Fe] 
and [Ca/Fe] is less than 0.03 dex in the thin disk. 

[Ti/Fe] was shown by EAGLNT to be a slowly decreas- 
ing function of [Fe/H] and the decrease continues to higher 
metallicity. Our data show a similar trend but the continu- 
ous decrease toward higher metallicity is less obvious with 
a comparatively large star-to-star scatter. There is no ev- 
idence that the scatter is correlated with Vlsr- We note 
that Feltzing & Gustafsson (1999) find a similar scatter in 
[Ti/Fe] for metal-rich stars with the Ti abundance based 
on 10-12 Til lines. 

7.2.3. Sodium and aluminum 



Na and Al are generally thought to be products of Ne and 
C burning in massive stars. The synthesis is controlled 
by the neutron flux which in turn depends on the ini- 
tial metallicity and primarily on the initial O abundance. 
Therefore, one expects a rapid increase of [Na/Mg] and 
[Al/Mg] with metallicity. But our data shows that both 
Na and Al are poorly correlated with Mg in agreement 
with EAGLNT. This means that the odd-even effect has 
been greatly reduced in the nucleosynthesis processes dur- 
ing the disk formation. 

When iron is taken as the reference element, we find 
that [Na/Fe] and [Al/Fe] are close to zero for [Fe/H] < 
-0.2, while EAGLNT found 0.1-0.2 dex differences be- 
tween [Fe/H] = -0.2 and [Fe/H] = -1 .0. Our results sup- 
port the old data by Wallerstein (1962) and Tomkin et al. 
( 1985D , who suggested [Na/Fe] ~ 0.0 for the whole metal- 
licity range of the disk stars. The situation is the same for 
Al; EAGLNT found an overabundance of [Al/Fe] ~ 0.2 
for [Fe/H] < —0.5, whereas we find a solar ratio for the 
low metallicity stars. As discussed in Sect. |5.3| this may, 
however, be due to a non-LTE effect. 

In the case of the more metal rich stars the abundance 
results for Na and Al are rather confusing. EAGLNT 
found that some metal-rich stars in the solar neighbour- 
hood are rich in Na, Mg and Al, but the existe nce o f 
such NaMgAl stars was rejected by Tomkin et al. ( 1997 ). 
Several further studies, however, confirmed the overabun- 
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Fig. 8. Abundance patterns for elements from O to Ba. The symbols are the same as in Fig. [?] and their size is 
proportional to stellar age. Note that the trends of [Al/Fe] and [K/Fe] may be spurious due to the neglect of non-LTE 
effects. 



dance of some elements again. Porte de Mortc (1996) 
found an overabundance of Mg but not of Na. Feltzing 
& Gustafsson (199E) confirmed the upturn of [Na/Fe] but 
their metal-rich stars did not show Mg and Al overabun- 
dances. In the present work we find a solar ratio of Na/Fe 



up to [Fe/H] ~ 0.1, and a rather steep upturn of [Al/Fc] 
beginning at [Fe/H] ~ —0.2. As discussed in Sect. our 
Al abundances may, however, be severely affected by non- 
LTE effects. We conclude that more accurate data on Na 
and Al abundances are needed. 
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7.2.4. Potassium 

[K/Fe] shows a decreasing trend with increasing metallic- 
ity for disk stars. The result supports the previous work 
by Gratton & Sneden ( 1987 ) but our data have a smaller 
scatter. Assuming that potassium is a product of explosive 
oxygen burning in massive stars, Samland (1998) repro- 
duces the observed trend rather well. Timmes et al. ( 1995| ) . 
on the other hand, predicts [K/Fe]< 0.0 for [Fe/H] < -0.6 
in sharp contrast to the observations. Given that the Kl 
resonance line at A7699, which are used to derive the 
K abundances, is affected by non-LTE as discussed in 
Sect. 5.3, it seems premature to attribute K to one of 



a elements. 

7.2.5. Vanadium, chromium and nickel 

V and Cr seem to follow Fe for the whole metallicity range 
with some star-to-star scatter. The scatter is not a result 
of mixing stars with different Vlsr, and the few very weak 
lines used to determine the abundances prevent us to in- 
vestigate the detailed dependence on metallicity and to 
decide if the scatter is cosmic or due to errors. 

Ni follows iron quite well at all metallicities with a 
star-to-star scatter less than 0.03 dex. Two features may 
be found after careful inspection. Firstly, there is a hint 
that [Ni/Fe] slightly decreases with increasing metallicity 
for —1.0 < [Fc/H] < —0.2. The trend is more clear, due to 
smaller star-to-star scatter, than in EAGLNT. Secondly, 
there is a subtle increase of [Ni/Fc] for [ Fc/H] > —0.2. 
Interestingly, Feltzing & Gustafsson (1999) found a slight 



increase of [Ni/Fe] towards even more metal-rich stars. 



7.2.6. Barium 

The abundance pattern of Ba is very similar to that of 
EAGLNT except for a systematic shift of about +0.07 dex 
in [Ba/Fe]. Both works indicate a complicated dependence 
of [Ba/Fe] on metallicity. First, [Ba/Fe] seems to increase 
slightly with metallicity for [Fe/H] < —0.7, and then keeps 
a constant small overabundance until [Fe/H] ~ —0.2, after 
which [Ba/Fe] decreases towards higher metallicities. 

Barium is thought to be synthesized by neutron cap- 
ture s-process in low mass AGB stars with an evolu- 
tionary timescale longer than that of iron-producing SNe 
la. Therefore, [Ba/Fe] is still slightly underabundant at 
[Fe/H] = —1.0. Ba is then enriched significantly at later 
stages of the disk evolution, but the decrease of [Ba/Fe] 
for more metal-rich stars beginning with [Fe/H] ~ —0.2 is 
unexpected. 

Given that the low [Ba/Fe] for some stars may be re- 
lated to their ages, the relation of [Ba/Fe] vs. [Fe/H] at 
different age ranges was investigated (see Fig. ^ . In agree- 
ment with EAGLNT, the run of [Ba/Fe] vs. [Fe/H] in old 
stars with logr 9 > 0.9 (—8 Gyr) and 0.7 < logr 9 < 0.9 
shows a flat distribution for [Fe/H] < —0.3 and a negative 
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Fig. 9. Different relations of [Ba/Fe] vs. [Fe/H] for stars 
with different age ranges: logrg < 0.7 (filled circles), 
logrg = 0.7— 0.9 (crosses) and logrg > 0.9 (open circles). 

slope for [Fe/H] > —0.3. All young stars with logrg < 0.7 
(~5 Gyr) have [Fe/H] > —0.3 and a clear decreasing trend 
of [Ba/Fe] with [Fe/H] is seen. In addition, there is a hint 
of higher [Ba/Fe] for younger stars both in the interval 
—0.7 < [Fc/H] < —0.3, where [Ba/Fe] is constant and in 
the interval [Fe/H] > —0.3, where [Ba/Fe] is decreasing. 
This is consistent with the formation of young stars at a 
later stage of the disk when long-lived AGB stars have en- 
hanced Ba in the ISM. The flat [Ba/Fe] for [Fe/H] < -0.3 
may be explained by the suggestion of EAGLNT that the 
synthesis of Ba in AGB stars is independent of metallic- 
ity, i.e. that Ba shows a primary behaviour during the 
evolution of the disk. But the age effect alone cannot ex- 
plain the underabundant [Ba/Fe] in metal-rich stars, be- 
cause [Ba/Fe] decreases with metallicity for all ages after 
[Fe/H] ~ —0.3. One reason could be that s-element syn- 
thesis occurs less frequently in metal-rich AGB stars pos- 
sibly because the high mass loss finishes their evolution 
earlier. 

8. Concluding remarks 

One of the interesting results of this study is that the 
oldest stars presently located in the solar neighbourhood 
have Vlsr < — 50kms~ 1 . Hence, they probably originate 
from the inner disk having i? m < 7 kpc. This is not co- 
incidentally found in our study. The EAGLNT sample 
contains about 20 such stars. As shown in both works, 
these stars are generally more metal-poor than other stars 
and they show a larger spread in [Fe/H] and [Ba/H] than 
in [a/H] (see Fig. |). According to EAGLNT, they have 
higher [a/Fe] than other disk stars at a metallicity about 
[Fe/H] = -0.7. 

Considering these different properties, we suggest that 
they do not belong to the thin disk. Firstly, they are older 
(10-18 Gyr) than other stars. Secondly, if they are thin 
disk stars, it is hard to understand why stars coming 
from both sides of the solar annulus have lower metal- 
licity than the local region. Thirdly, these stars show a 
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relatively small metallicity dispersion at such early Galac- 
tic time, i.e. smaller than stars at 8-10 Gyr. This is not in 
agreement with the effect of orbital diffusion working dur- 
ing the evolution of the thin disk, which suggests larger 
metallicity dispersion for older stars. Finally, the Wlsr 
dispersion of these stars is about 40 kms" 1 , considerably 
larger than the typical value of about 20 kms" 1 for thin 
disk stars. Consistently, the kinematics, age, metallicity 
and abundance ratios of these stars follow the features of 
the thick disk: Vlsr < -50 kms" 1 , ct(WlSr) ^ 40 kms" 1 , 
r > 10 Gyr, [Fe/H] < -0.5 and [a/Fe] ~ 0.2. We conclude 
that these oldest stars in both EAGLNT and this work 
are thick disk stars. Hence, they are probably not result- 
ing from an inside-out formation of the Galactic disk, but 
have been formed in connection with a merger of satellite 
components with the Galaxy. 

Concerning the abundance connection of the thick disk 
with the thin disk, our data for [a/Fe], shown in Fig. |io| , 
suggest a more smooth trend than those of EAGLNT, who 
found a correlation between [a/Fe] and R m at [Fe/H] ~ 
—0.7. We leave the issue open considering the small num- 
ber of these stars in our work. Two stars marked by their 
names in Fig. [l^ may be particularly interesting because 
they show significantly higher [a/Fe] than other stars. 
Fuhrmann & Bernkopf (199£) suggest that one of them, 
HD 106516A, is a thick-disk field blue straggler. It is un- 
clear if this can explain the higher [a/Fe]. HD 97916 is 
a nitrogen rich binary (Beveridge & Sneden 1994) with 
C^lsr = — 117kms~ 1 and TVlsr = lOlkms" 1 (typical for 
halo stars), but with Vlsr = 22 kms" 1 similar to the value 
for thin disk stars. Surprisingly, this star is also very young 
(5.5 Gyr) for it's metallicity. 

It is interesting to re-inspect the observational results 
for thin disk stars excluding the thick disk stars. More di- 
rect information on the evolution of the Galactic thin disk 
will be then obtained. In summary, the thin disk is younger 
(not older than 12 Gyr), more metal-rich ([Fe/H] > —0.8) 
and has a smaller [a/Fe] spread (0.1 dex) without the 



mixture of the thick disk stars. In particular, the AMR is 
more weak and there seems to exist a radial metallicity 
gradient. All these features agree better with the present 
evolutionary models for the Galactic disk. 

We emphasize here that there is no obvious gradient 
in [a/Fe] for the thin disk at a given metallicity. Such 
a gradient was suggested by EAGLNT based on higher 
[a/Fe] of the oldest stars with R m < 7 kpc than stars with 
R m > 7 kpc (see their Fig. 21). After we have ascribed 
these oldest stars to the thick disk, the abundance gradient 
disappears. 

Our study of relative abundance ratios as a function of 
[Fe/H] suggests that there are subtle differences of origin 
and enrichment history both within the group of a ele- 
ments and the iron-peak elements. Nucleosynthesis the- 
ory predicts that Si and Ca are partly synthesized in SNe 
la, while O and M g are only produced in SNe II (Tsuji- 
moto et al. 1995). Our data suggest, however, that SNe 



la may also be a significant synthesis site of Mg, because 
[Mg/Fe] shows a trend more similar to [Si/Fe] and [Ca/Fe] 
than to [O /Fe] . Ti may not lie in a smooth extension of Si 
and Ca, because there is a hint of a decrease of [Ti/Fe] 
for [Fe/H] > —0.4 not seen in the case of Si and Ca. 
The situation for the odd-Z elements is more complicated. 
The available data for Na and Al show confusing disagree- 
ments; EAGLNT finds an overabundance of 0.1 to 0.2 dex 
for [Na/Fe] and [Al/Fe] among the metal-poor disk stars, 
whereas our study points at solar ratios. Two other odd-Z 



elements, K and Sc (Nissen et al. 1999), behave like a ele 



ments, but the result for K is sensitive to the assumption 
of LTE. The iron-peak elements also show different be- 
haviours: V, Cr and Ni follow Fe very well, while [Mn/Fe] 
(Nissen et al. 199E ) decreases with decreasing metallicity 
from [Mn/Fe] ~ 0.0 at [Fe/H] = 0.0 to [Mn/Fe] ~ -0.4 at 
[Fe/H] = —1.0. We conclude that the terms "a elements" 
and "iron-peak elements" do not indicate productions in 
single processes, and that each element seems to have a 
unique enrichment history. 
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